第 27 卷 第 2 期 工 程 数 X E 报 Vol. 27 No. 2 
2010404 H CHINESE JOURNAL OF ENGINEERING MATHEMATICS Apr. 2010 


文章 编号 :1005-3085(2010)02-0313-08 
马 氏 链 预测 方法 的 统计 试验 研究 * 


夏 乐 天 !， UE 
(1- 河 海 大 学 理学 院 ， 南 京 210098; 2- 河 海 大 学 水 资源 环境 学 院 ， 南 京 210098) 

搞 要 : 近年 米 ， 基 于 绝对 分 布 的 马尔 可 夫 链 预测 (ADMCP) 方 法 、 琶 加 马尔 可 夫 链 预测 (SPMCP) 方 
法 和 加 权 马 尔 可 夫 链 预测 (WMCP) 方 法 在 各 种 水 文 序 列 的 预测 理论 中 得 到 了 广泛 的 应 用 。 
尽管 各 种 马尔 可 夫 链 预测 方法 在 有 关 应 用 实例 中 效果 都 不 错 ， 但 是 ， 实 际 上 它们 的 预测 效果 
是 有 差别 的 。 本 文通 过 统计 试验 手段 对 这 三 种 马尔 可 夫 链 预测 方法 进行 了 比较 研究 ， 结 论 表 
Jj: WMCP 方法 预测 精度 最 高 ，SPMCP 方法 预测 精度 次 之 ，ADMCP 方法 预测 精度 相对 来 讲 


最 差 。 
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1 引言 


近年 来 ， 张 宕 等 叫 研究 了 马尔 可 夫 链 预测 理论 在 矿区 降水 灾害 预测 中 的 应 用 ; 郑 文 瑞 
FA 研究 了 马尔 可 夫 链 预测 理论 在 水 污染 状态 风险 评价 中 的 应 用 ， 汉 利 华 等 由 用 马尔 可 夫 链 预 
测 理论 研究 了 区 域 干 旱 的 变化 趋势 ， 刘 德 辅 等 向 用 马尔 可 夫 链 预测 理论 研究 了 沿海 和 河口 城市 
防 灾 设 防 标准 的 系统 分 析 :， 宋 印 胜 回 和 王 渺 林 [ 研究 了 马尔 可 夫 链 理论 在 水 位 预测 中 的 应 用 ; 
潘 俊 等 中 研究 了 马尔 可 夫 链 理论 在 水 库 调节 运行 中 的 应 用 ， 钱 家 忠 等 回 在 地 下 水 资源 评价 降水 
量 的 分 析 中 运用 了 马尔 可 夫 链 预测 理论 ; 文献 [9] 研究 了 地 区 洪水 的 马尔 可 夫 链 预测 ， 汉 耀 龙 
等 上 9 研究 了 马尔 可 夫 链 在 河流 丰 枯 状况 预测 中 的 应 用 ， 孙 才 志 等 上 研究 了 马尔 可 大 链 在 降水 
丰 枯 状况 预测 中 的 应 用 ;， 夏 乐天 等 外 研究 了 马尔 可 夫 链 在 太湖 流域 和 长 江 中 下 游 地 区 的 梅 南 
预报 中 的 应 用 等 。 在 所 发 表 的 大 量 研究 论文 中 ， 各 种 马尔 可 夫 链 预测 Da 方法 都 有 不 同 程度 的 
应 用 ， 但 是 常用 的 三 种 马尔 可 夫 链 预测 方法 的 优 劣 性 却 一 直 鲜 有 讨论 。 

由 于 我 国有 关 水 文 计 算 的 规范 规定 : 我 国 主要 河流 的 径流 量 和 降水 量 均 可 假定 为 服 
从 P-II 型 分 布 或 三 参数 对 数 正 态 分 布 。 因 此 本 文 的 主要 目的 是 利用 随机 模拟 技术 ， 生 成 服 
从 P-II 型 分 布 和 三 参数 对 数 正 态 分 布 的 相依 伪 随 机 数 序列 ， 然 后 再 对 各 种 不 同 的 马尔 可 夫 链 
预测 方法 进行 比较 分 析 并 讨论 各 种 因素 对 预测 精度 的 影响 。 所 讨论 的 主要 问题 是 : 

1) 各 种 马尔 可 夫 链 预测 方法 的 适应 性 ， 包 括 对 不 同 总 体 分 布 的 适应 性 (变量 总 体 分 布 为 对 
数 正 态 与 P-II 型 )， 和 对 不 同 的 序列 生成 模型 的 适应 性 ; 

2) 讨论 各 种 马尔 可 夫 链 预测 方法 的 有 效 性 ， 得 出 哪 种 预测 方法 较 优 的 结论 。 


2 试验 设计 
2.1 伪 随 机 数 序列 的 长 度 和 分 布 参数 的 设计 
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由 于 天 然 河流 的 实测 水 文 资料 比较 短 ， 一 般 少 于 60 年 ， 因 此 只 有 在 小 样本 时 讨论 预测 方法 
的 优 劣 ， 才 会 有 实际 的 水 文 意义 。 出 于 这 种 考虑 ， 本 文 对 于 生成 的 序列 ， 仅 考虑 长 度 为 100 的 
服从 P-II 型 分 布 和 三 参数 对 数 正 态 分 布 的 相依 伪 随 机 数 序列 ， 利 用 适当 的 分 组 方法 确定 这 些 
数据 的 状态 。 假 定 实测 资料 的 长 度 为 60 年 ， 利 用 这 60 年 的 资料 所 对 应 的 状态 预测 出 未 来 40 年 
的 状态 ， 然 后 与 这 40 年 的 理论 状态 做 比较 ， 为 消除 偶然 性 ， 可 由 生成 的 序列 ， 连 续 考虑 长 度 
为 100 的 服从 P-III 型 分 布 和 三 参数 对 数 正 态 分 布 的 相依 伪 随 机 数 序列 ， 分 别 利 用 前 60 年 的 资 
料 所 对 应 的 状态 预测 出 未 来 40 年 的 状态 ， 然 后 与 后 40 年 的 理论 状态 做 比较 ， 算 出 平均 误差 ， 
最 后 定 出 各 种 马尔 可 夫 链 预测 方法 的 优 劣 。 为 在 一 定 的 参数 范围 内 讨论 评价 各 种 马尔 可 夫 链 预 
测 方法 的 优 劣 ， 共 采用 了 5 组 总 体 参 数 ， 见 表 1。 


表 1: 总 体 参数 值 表 


经 验 表明 : 大 多 数 实 际 径流 量 和 降水 量 资料 的 参数 是 在 表 1 的 参数 值 的 范围 内 ， 序 列 的 一 
阶 相关 系数 统一 规定 为 0.2， 由 于 一 般 的 水 文 序列 均 为 弱 相 关 序 列 。 

2.2 平稳 马尔 可 夫 链 的 模拟 

许多 水 文 过程 都 表现 出 明显 的 序列 相关 ， 通 常 以 ps(k) 表示 随机 变量 XX 的 上 阶 序列 相关 系 
数 ， 如 果 可 以 把 ps(k) 近似 地 表示 成 (pz(1)*， 则 随机 变量 匀 的 时 间 序 列 可 以 用 一 阶 平稳 马尔 
可 夫 链 作为 模型 4。 

1961 年 ， 在 径流 调节 的 研究 中 开始 使 用 一 阶 Markov 模 型 ， 即 考虑 序列 相 邻 项 的 相关 关系 
以 生成 水 文 序列 中。 直到 目前 ， 这 种 模型 仍 是 水 文学 中 使 用 最 多 的 一 种 ， 故 本 节 对 该 种 模型 
生成 的 序列 给 出 较为 详尽 的 介绍 。 

确切 地 说 ， 一 阶 马尔 可 夫 链 由 下 式 定义 


Zi+1 = Hz + pz(1) (Ti — uz) + Eiti, (1) 
其 中 z; 是 马尔 可 夫 链 在 时 刻 i 处 的 值 , ys 是 XX 的 均值 , ps(1) 是 一 阶 序列 的 自 相关 系数 , (6i) 是 
满足 Ec; = 0, Vare; = o2 的 独立 随机 变量 序列 。 这 一 模型 表明 : 义 在 某 一 时 刻 的 取 值 ， 只 与 
其 前 一 时 刻 的 取 值 有 关 ， 再 加 上 一 个 随机 变量 。 设 以 o2 表示 义 的 方差 ， 考 虑 到 ei41 与 X; 相 

互 独立 ， 故 可 得 2 与 02 的 关系 如 下 

oz = Var(Xit1) = E[Xi — us]? = 920) E(X; — uz)? + Ee, = p2(1)o2 + 02, 
即 

ez =02[1— pz(1)]. (2) 
假设 无 的 分 布 是 正 态 的 NUue,a23)， 则 由 (1) 可知 的 分 布 是 N(0,02)。 设 以 + 表示 标 准 
正 态 随 机 变量 ， 则 ros V1 一 P2(1) 是 服从 分 布 N(0,o2) 的 随机 变量 。 于 是 要 抽样 服从 分 
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布 W(pz;ca2)、 一 阶 序列 的 自 相 关系 数 是 pe(1)、 且 符合 一 阶 Markov 模 型 的 序列 {Xi}， 其 抽 
样 公式 应 如 下 式 所 示 

Titl = Hs + pz(l)(Ti— Hz)+ritioz V1 — p2(1). (3) 


为 利用 上 式 生 成 Til» 首先 要 使 用 某 种 方法 得 出 Hx: Ox, px (1) 的 估计 和 值 x, $5; rx (1). i, Sr 
的 估计 值 皆 为 已 知 ，rz(1) 可 按 自 相关 系数 的 定义 采用 下 式 计算 


n—1 n—1 n—1 


zu Dr TiTit+l 一 Guy 2 Ti x Til 
i i i=1 


rl) = 一 二 z 页， 
[4 Esto eim (o) ] [tnm 17 
bi TiTitl 一 TPE z x 
- AL o T m (4) 
[Esto OR YT Dm S9 TA 


和 的 初 值 zz， 可 以 由 分 布 NW(uc;,c3s) 抽 得 。 为 了 消除 初 值 的 影响 ， 可 以 将 生成 序列 的 前 几 
项 (例如 前 50 项 至 100300) 去 掉 。 

式 (3) 曾 被 广泛 地 用 于 生成 年 径流 序列 。 但 由 于 正 态 分 布 可 取 负 值 的 特性 与 径流 过 程 的 特 
性 (不 可 能 出 现 负 值 ) 有 差异 ， 特 别 是 us 较 小 而 os 较 大 时 ， 用 (3) 式 生成 的 zi 常会 出 现 负 值 。 
如 负 值 出 现 的 机 会 不 多 ， 则 可 只 用 这 负 的 x; 来 生成 下 一 个 zir1， 然 后 将 xz; 从 序列 中 去 掉 舍弃 
不 用 。 但 当 负 的 zi 出 现 很 多 时 ， 则 表明 径流 过 程 X, 的 分 布 ， 与 正 态 分 布 差异 显著 ， 此 时 必须 
考虑 其 它 分 布 。 在 水 文 统 计 中 常用 的 分 布 有 两 种 ， 即 Gamma 分 布 或 P-II 型 分 布 及 对 数 正 态 


WEERA X 能 服从 P-II 型 分 布 时 ， 可 使 用 下 式 1g 
Tit1 = Ha + pa(1)(Ti — Hz) + i102 V/A — p2(1), (5) 

其 中 : 
saca (re ER LR (6) 


Tit1 N(0,1) BREL Cos 为 € 的 偏 态 系数 ， 由 下 式 给 出 


1-p2(1) 
Cs; = ir-20g (7) 


为 使 生成 的 天 能 服从 三 参数 对 数 正 态 分 布 ， 可 假定 Y = In(X —a) 服从 正 态 分 布 (uy, 02), 
于 是 Y 可 采用 (3) 式 的 生成 方法 ， 即 


Yiti = Hy + py (1)(yi — ny) + Ti+10y V1 — p2(1), (8) 
HEP uy, oy, py(1) 分 别 为 Y 的 均值 、 均 方差 和 一 阶 序 列 的 自 相关 系数 。 直 接 用 z 表示， 这 一 生 
成 方法 为 
zi = Qa+{exp[py(l — py(1))]} Gri — oa) O61, (9) 
其 中 
6441 = exp ([1 — p2 (1)] oyri }- 
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具体 计算 时 ， 首 先 用 某 种 方法 估计 出 瑟 的 有 关 参 数 Hr, oz, Cs,, pz(1) GEP Cos XX KI 
态 系数 )， 再 由 天 的 前 二 阶 矩 与 参数 间 的 关系 19， 可 得 如 下 几 个 公式 


Oy = vV/In(1 一 72), a = Lia (1 T C, /1]), 
[zur cuo THEN T 
2 2 ? 


7 = 
1 2 1 g? 
Hy =ln(ps —a)- 51-90) py(1) = zz In[(e?v — 1)pz(1) * 1]. 
y 


将 这 些 算得 的 参数 代入 (8) SX EX (9) I0, RU Pr ERKI ZU (Xi), ERE IR RE X H a E S 
数 : Hax, Ox, Cs, 以 及 pz(1)。 


3 ”计算 方案 的 设计 与 统计 试验 成 果 的 综合 分 析 


31 计算 方案 、 预 测 偏差 分 析 与 统计 试验 成 果 

3.1.1 计算 方案 与 预测 偏差 

根据 (5) 式 和 (9) 式 ， 按 表 1 中 给 定 的 各 组 参数 分 别 生 成 长 度 为 100000 的 服从 P-II 型 
分 布 和 对 数 正 态 分 布 的 马尔 可 夫 链 序列 mm ， 依 次 取 100 个 数据 (分 别 考虑 取 50 次 、200 次 
和 1000 次 )， 用 三 种 马尔 可 夫 链 预测 方法 分 别 去 做 预测 。 即 假定 实测 资料 的 长 度 为 60 年， 用 
文献 [13] 3.3.1 中 介绍 的 分 级 方法 ， 对 一 组 给 定 的 参数 ai, az, as 和 as， 把 这 60 年 的 指标 值 分 
级 ， 利 用 这 60 年 的 资料 所 对 应 的 状态 分 别 用 三 种 马尔 可 夫 链 预测 方法 预测 出 未 来 4 年 的 状 
态 ( 共 做 50 次 、200 次 和 1000 次 );， 然后 把 逐步 预测 出 的 状态 与 理论 状态 作 差 ， 对 这 个 差 的 绝 
对 值 先 按 次 数 作 均 值 和 均 方 差 ， 代 表 逐 步 预 测 偏差 和 逐步 预测 偏差 均 方 养 ， 再 考虑 逐步 预测 偏 
差 和 逐步 预测 偏差 均 方 差 的 平均 值 ， 可 称 为 预测 偏差 均值 和 预测 偏差 均 方 差 ， 这 个 预测 偏差 均 
值 和 预测 偏差 均 方 差 就 是 度量 预测 方法 优 劣 的 两 个 标准 。 显 然 ， 预 测 偏差 均值 和 预测 偏差 均 方 
差 的 值 是 越 小 越 好 ， 它 们 反映 了 预测 的 精度 。 这 个 过 程 可 用 数学 语言 表述 如 下 。 

对 于 100 个 生成 的 服从 某 种 分 布 (P-THI 型 分 布 或 对 数 正 态 分 布 ) 的 数据 zi,，…. ,z6o,…… ,zioo， 
把 后 面 40 个 数据 对 应 的 状态 记 为 z1,… , zao， 作 为 理论 数据 , 利用 前 面 60 个 数据 z1,.… ,zeo 所 
对 应 的 状态 ， 通 过 某 种 马尔 可 夫 链 预测 方法 逐步 预测 出 后 面 40 个 数据 所 对 应 的 状态 记 
为 2 四 ,… ,z 久 ,7 二 1,… ,c， 其 中 为 模拟 次 数 ， 则 第 i 步 逐 步 预 测 偏差 定义 为 


msl- i=1,...,40, (10) 
预测 偏差 均值 定义 为 
_ 1 40 z 1 40 c G) 
Wr -= žij, (11) 
1 ~ j B 
ao d EP som) , 1 一 1 ,40, (12) 
预测 偏差 均 方 差 定义 为 


12 
8 — 40 2 : (13) 
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3.1.2. 预测 成 果 汇 编 


取 定 EX = 1000，pz(1) = 0.2, (01,02,03, 04) = (1.2,0.6,0.4,1.0)， 把 两 种 分 布 的 预测 统计 
试验 简要 结果 整理 如 下 ， 见 表 2 和 表 3。 


表 2: 了 P-III 型 分 布 预测 统计 试验 简要 乡 


5 
%7 


模拟 次 数 


0. 

0.8190 
0.8188 
0.7980 


LEGO | ow | oo | as | 215 | s | 


3.2. ”统计 试验 成 果 分 析 

3.2.1 三 种 马尔 可 夫 链 预测 方法 的 预测 精度 不 同 

就 参数 (C,, C.) 的 三 组 值 的 平均 情况 来 说 ， 三 种 马尔 可 夫 链 预测 方法 的 预测 精度 是 不 同 
的 ( 见 表 2、 表 3 中 各 部 分 的 最 后 二 行 )， 具 体 地 说 : 加 权 马 尔 可 夫 链 预测 方法 的 预测 精度 最 
高 ， 和 盈 加 马尔 可 夫 链 预测 方法 的 预测 精度 次 之 ; 基于 绝对 分 布 的 马尔 可 夫 链 预测 方法 的 预测 精 
度 最 低 。 i 

事实 上 ， 对 P-II 型 分 布 而 言 ， 以 基于 绝对 分 布 的 马尔 可 夫 链 预测 方法 (ADMCP 法 ) 为 基 
准 ， 肥 加 马尔 可 夫 链 预测 方法 (SPMCP 法 ) 的 预测 偏差 均值 要 提高 4.78 ~ 4.81 个 百分点 ， 预 测 
偏差 均 方差 要 提高 1.54 ~ 2.15 个 百分点 ， 加 权 马 尔 可 夫 链 预测 方法 (WMCP 法 ) 的 预测 偏差 均 
值 要 提高 5.80 ~ 6.09 个 百分点 ， 预 测 偏差 均 方 差 要 提高 2.53 ~ 3.00 个 百分点 。 
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对 三 参数 对 数 正 态 分 布 而 言 ， 以 基于 绝对 分 布 的 马尔 可 夫 链 预测 方法 (ADMCP 法 ) 为 
基准 ， 重 复试 验 200 次 以 上 ， 炙 加 马尔 可 夫 链 预测 方法 (SPMCP 法 ) 的 预测 偏差 均值 要 提 
高 4.69 ~ 5.54 个 百分点 ， 预 测 偏差 均 方 差 要 提高 1.71 ~ 2.19 个 百分点 ; 加 权 马 尔 可 夫 链 
预测 方法 (WMCP 法 ) 的 预测 偏差 均值 要 提高 5.82 ~ 6.92 个 百分点 ， 预 测 偏 差 均 方差 要 提 
高 3.14 ~ 3.51 个 百分点 。 


表 3: 对 数 正 态 分 布 预测 统计 试验 简要 结果 


s| = 


1.50 
2.00 
2.50 


1000 
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Coson | 


MaRa om | 000 | se | 


3.2.2 ”统计 试验 的 次 数 c 对 预测 结果 有 一 定 的 影响 

事实 上 试验 的 次 数 过 少 ， 例 如 c = 50 时， 预测 精度 较 低 ，c > 200 时 ， 预 测 精度 趋向 稳定 。 

3.2.3 ”各 种 马尔 可 夫 链 预测 方法 的 预测 精度 具有 对 不 同 分 布 的 适应 性 

以 试验 次 数 c = 200 为 例 ， 基 于 了 P-II 型 分 布 的 五 组 参数 的 平均 预测 偏差 均值 为 ADMCP 法 
0.8401，SPMCP 法 0.7997，WMCP 法 0.7914, 基于 对 数 正 态 分 布 的 五 组 参数 的 平均 预测 偏差 均 
值 为 ADMCP 法 0.8526，SPMCP 法 0.8054，WMCP iX: 0.7936. 

基于 P-IIJ 型 分 布 的 五 组 参数 的 平均 预测 偏差 均 方 差 为 ADMCP 法 0.8132, SPMCP i; 0.7957, 
WMCP 法 0.7888， 基 于 对 数 正 态 分 布 的 五 组 参数 的 平均 预测 偏差 均 方差 为 ADMCP 法 0.8162， 
SPMCP 法 0.7983，WMCP 法 0.7875。 可 见 基 于 两 种 分 布 的 预测 精度 没有 本 质 的 差别 。 
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3.2.4 参数 (Cu, C) 的 不 同 对 预测 精度 的 影响 显著 

一 般 来 说 ， 随 着 C, 和 Cs 的 增 大 ， 预 测 偏差 呈 下 降 趋势 ， 且 WMCP 法 的 预测 结果 最 为 稳 
健 。 

以 试验 次 数 c = 200 为 例 ， 基 于 了 -IJ 型 分 布 的 五 组 参数 ，ADMCP 法 的 max mz = 0.9487, 
min mz = 0.7322, max sz = 0.8372, min sZ = 0.7945; SPMCP 法 的 max mz = 0.8975, min 
mz = 0.6849, max sz = 0.8255, min sz = 0.7672; WMCP 法 的 max mz = 0.8766, min mz = 
0.6898, max sz = 0.8085, min sz = 0.7616. mz TR Æ 1E 0.1868 ~ 0.2165, sz 的 极 差 在 0.0469 
~ 0.0583， 一 般 说 来 ，ADMCP 法 的 极 差 大 一 些 ，SPMCP 法 的 极 差 略 小 ，WMCP 法 的 的 极 差 
最 小 ， 这 说 明 : 在 参数 (Cu,C。) 变化 时 ，WMCP 法 的 预测 结果 最 为 稳健 。 

基于 对 数 正 态 分 布 的 五 组 参数 , ADMCP 法 的 max mz = 0.9315, min mz = 0.7462, max sZ 
— 0.8379, min sz — 0.7878; SPMCP 法 的 max mz = 0.8693, min mz — 0.7135, max sz = 
0.8135, min sz = 0.7763, WMCP 法 的 max mz = 0.8568, min mz = 0.6992, max sz = 
0.8042, min sz = 0.7658. mZ AJR 2: 1E 0.1558 ~ 0.1863, sz HJ TR2E1E0.0384 ~ 0.0501, — 
般 说 米 也 是 ，ADMCP 法 的 极 差 大 一 些 ，SPMCP 法 的 极 差 略 小 ，WMCP 法 的 极 差 最 小 ， 这 说 
8j: 在 参数 (Cu, C,) 变化 时 ，WMCP 法 的 预测 结果 也 最 为 稳健 。 

由 于 篇 幅 所 限 ， 我 们 将 另外 撰文 进一步 讨论 Cu 固定 ，C 逐步 增加 或 Cs 固定 ，C 逐步 增加 
对 预测 精度 的 影响 ， 以 及 样本 容量 、 一 步 自 相关 系数 pa(1)、 指 标 值 分 级 的 不 同 对 预测 精度 的 
影响 。 


4 小结 


本 文通 过 统计 试验 对 各 种 不 同 的 马尔 可 夫 链 预测 方法 进行 了 比较 分 析 并 讨论 了 各 种 因素 对 
预测 精度 的 影响 。 

统计 试验 结果 表明 : 

1) 加 权 马 尔 可 夫 链 预测 方法 的 预测 精度 最 高 ， 受 加 马尔 可 夫 链 预测 方法 的 预测 精度 次 
之 ; 基于 绝对 分 布 的 马尔 可 夫 链 预测 方法 的 预测 精度 最 低 。 

2) 统计 试验 的 次 数 c 对 预测 结果 有 一 定 的 影响 ， 当 c > 200 时 ， 预 测 精度 趋向 稳定 。 

3) 各 种 马尔 可 夫 链 预测 方法 的 预测 精度 具有 对 不 同 分 布 的 适应 性 。 

进一步 讨论 可 参见 文献 [13]。 
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Study on Statistical Experiment of Markov Chain 
Prediction Methods 


XIA Le-tian!, ZHU Yuan-shen? 


(1- College of Science, Hohai University, Nanjing 210098; 

2- College of Hydrology and Water Resources, Hohai University, Nanjing 210098) 
Abstract: In recent years, the weighted Markov chain prediction method (WMCP), the Markov chain 
prediction method based on absolute distribution (ADMCP) and probability summation (SPMCP) 
have been widely used in the various hydrology series. Although the application results of various 
Markov chain prediction methods in real situations are seemed good, the prediction accuracy of these 
three kinds of Markov chain prediction methods is different. The three prediction methods are studied 
by using statistical experiment. The experiment results indicate that the prediction accuracy of the 
WMOP method is the highest, the ADMCP method is the lowest, and the SPMCP method is in the 
middle. 
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